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ABSTRACT 


Observations by the Fermi Large Area Telescope of 7 -ray millisecond pulsar 
light curves imply copious pair production in their magnetospheres, and not ex¬ 
clusively in those of younger pulsars. Such pair cascades may be a primary source 
of Galactic electrons and positrons, contributing to the observed enhancement 
in positron flux above ~10 GeV. Fermi has also uncovered many new millisec¬ 
ond pulsars, impacting Galactic stellar population models. We investigate the 
contribution of Galactic millisecond pulsars to the flux of terrestrial cosmic-ray 
electrons and positrons. Our population synthesis code predicts the source prop¬ 
erties of present-day millisecond pulsars. We simulate their pair spectra invoking 
an offset-dipole magnetic field. We also consider positrons and electrons that 
have been further accelerated to energies of several TeV by strong intrabinary 
shocks in black widow and redback systems. Since millisecond pulsars are not 
surrounded by pulsar wind nebulae or supernova shells, we assume that the pairs 
freely escape and undergo losses only in the intergalactic medium. We compute 
the transported pair spectra at Earth, following their diffusion and energy loss 
through the Galaxy. The predicted particle flux increases for non-zero offsets of 
the magnetic polar caps. Pair cascades from the magnetospheres of millisecond 
pulsars are only modest contributors around a few tens of GeV to the lepton 
fluxes measured by AMS—02, PAMELA, and Fermi, after which this component 
cuts off. The contribution by black widows and redbacks may, however, reach 
levels of a few tens of percent at tens of TeV, depending on model parameters. 

Subject headings: cosmic rays — pulsars; general — stars: neutron 
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Introduction 


Telescope (LAT; 


Aguilar et al 


Ackerman u et al 


2013 


2014 : 


2012 


Accardo et al 


Recent mea surements by PAMELA flAdriani et al.ll2009l. 120131) . Fermi Large Area 


, and the Alpha Magnetic Spectrometer {AMS—02; 


201411 have provided hrm evidence that the positron 


fraction (PF) (^{e'^)/[(j){e'^) + 0(e“)], with 0 the flux, is increasing with energy above 
~ 10 GeV. Improved spectral measurements for 30 months of AMS—02 data extended the 
PF up to 500 GeV, and indicated a leveling off of this fraction with energy, as well as the 
PF being consistent with isotropy. 

Secondary positrons are created during inelastic collisions between cosmic-ray nuclei 
and intergalactic hydrogen, which produce charged pions that in turn decay into positrons, 
electrons, and neutrinos. The fraction of this secondary component with respect to the 
total (electron -|- positron) cosmic-ray spectrum is expected to smo othly decrease with 


Moskalenko fc Strong 


energ y within the standard framework of cosmic-ray transport (e.g., 

19981 ) However, the AMS— 02 electron spectrum is softer than the positron one in the 


range 20 — 200 GeV flAguilar et al 


2014 1. and the measured PF rises with energy, pointing 
to nearby sources of primary positron^. Moreover, the rising PF can be ascribed to a 


^This, however, depends on model assumptions, i.e., a concave electron spectrum may 
lead to a rising PF with energy. 


^Such an additional source o 
lation origin (e.g., 
or of astrophysica 


Grasso et al. 


primary positrons may be either of da r k matter annih i- 


2009 


Fan et al. 


20101 : 


Porter et al. 


micro quasar jets fl Gupta fc Torres 


origin, including su 


sar wind nebulae (e.g., 


3ernovae (e.g.. 


20 


Blasi fc Amato 


Blasi 


rnolecular cloud s (iDogiel fc Sharov 


2009 


2011 


Lin et al 


Delahave et al. 


2015 1. 


201011 . 


2011 


Serpico 


Arons 

1981; 

Harding V 

Ramatv 

1987: 

Boulares 

1989; 

Ghi et al. 

1996; 

Zhang V Gheng 

2001 

Grimani 

2 OO 7 I 

ProfumJ 

2 OI 2 I 

Kistler & Yiiksel 

2009; 

Hooper et al. 

2009; 

Yiiksel et al. 

2009 


1990 1. pul- 


201211. young or mature pulsars (e.g. 
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hardening of the positron spectrum (up to 200 GeV, after which it softens with energy), 
and not a softening in electron spectrum above 10 GeV. 

Alternatively, it has been argued that the observed rise in PF with energy may be 


explained purely by secondary positrons originating in the interste 


without the need to invoke a primary positron source. 


Shaviv et al 


l ar me dium (ISM), 
(1200911 demonstrated 


that an inhomogeneous distribution of supernova remnants (SNRs), such a s a strong 


conce ntration in the Galactic spiral arms, may explain the PF shape (see also 


Gaggero et al. 


2014J . who note that an unrealistically steep index for the primary electron spectrum needs 


to be invoked when assuming a homogeneous or smoothly varying source distribution; 
however, they do find eviden ce for an ext r a / se condary charge-symmetric electron-positron 
source to explain the data). iMoskalenkol (120131) p ointed out that the concave shape of the 


primary 
the PF. 


electron spectrum of 


Shaviv et al 


(1200911 introduces an arguably artihcial rise in 


Gowsik fc Burchl (120101) put forward a model assuming that a signihcant fraction 


of the boron below 10 GeV is generated through spallation of cosmic-ray nuclei in small 
regions around the sources. In this case, the contribution from spallation in the ISM would 


have a flat or weak energy dependence, and the GeV posi 
be generated through cosmic-ray interactions in the ISM. 


rons would almos t exclusively 


Moskalenko ( 2013 1 noted that 


such sources should be observable as very bright GeV y-ray sources with soft spectra, while 
the diffuse emission would be signihcantly dimmer than observed. This sce nario is also at 


odds with current estimates of the supernova birth rate. 


Rlnm et al 


(120131) found an upper 


bound to the positron flux by neglecting energy losses, arguing that the flattening of the 
PF seen by AMS-02 around several hundred GeV is consistent with a purely secondary 


Gendelev et al. 

2010 ; 

Yin et al. 

2013 


Feng fc Zhang 


formed 

(MSPs; 


3 V the merger of two white dwarfs (iKashivama et al. 


201511, ‘whi te dwarf pulsars’ mainly 


Kisaka fc Kawanaka 


201 ll) . and millisecond pulsars 


201211 . In this paper, we investigate the latter source class. 
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origin for the positrons. iMoskalenkol (120131) noted that their arguments imply quite hard 


injection spectra for primary nuclei, in contradiction to 7 -ray observations of SNRs that 
seem to require rather steep spectra. In addition, a very fast escape time for the positrons 



requiring a primary source of positrons in this case. 

MSPs are the oldest population of rotation-powered pulsars, characterized by low 


surface magnetic helds, and are thought to hav e acquired t 


spin-up by accretion from a binary companion flAlpar et ah 


reir ve ry short periods through 


19821) . For the most part, they 


have not been considered as an important source of cosmic-ray positrons since the majority 


lie below the c 


Harding et ah 


eath l i nes for high-multiplici ty pair cascades (assuming dipole magnetic helds; 


(iHarding et ah 


2002 


Zhang &: Chen 


20031) and were thus considered to be pair-starved 


20051 ). However, this picture cha nged with t 


le de 


emission from a large number of MSPs by Fermi (lAbdo et al 


section of pulsed 7 -ray 


2013 ). Most of the 7 -ray light 


curves show narrow, double peaks trailing the radio peaks, very similar to tho se of younger 


pulsars. Such light curves can only be ht by outer magnetospheric gap models (jVenter et al. 


2009 


JphnsqnetaL 


2 OI 4 J ). The existence of narrow accelerator gaps requires large numbers 


of electron-positron pairs (high multiplicity) to screen the electric held parallel to the 
magnetic held in the open magnetosphere interior to (at lower colatitudes than) the gaps. 


It has been suggested that distortions of the surface magnetic hel d may increase 


pair 


production for MSPs, eit her in the form of higher multipoles (e.g. 


or ohse t polar caps (PCs; 


Arons 


1996 


Harding fc Muslimov 


2 nila 


Zhang fc Cheng 


M). 


29931). 


Harding fc Muslimov 


fj2911al ) found that even small ohsets of the PC from the magnetic axis (a small fraction of 
the stellar radius) can greatly enhance the pair multiplicity. This is due to the increase in 
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accelerating electric field on one side of the PC, which stems from the decrease in curvature 
radius of the distorted magnetic held. Furthermore, MSPs produce pairs with energies 
around 100 times higher that those of young pulsars, due to their relatively low magnetic 
helds which require a higher photon energy for magnetic p hoton pair production to tak e 


place. In this case, the pair spectra extend to several TeV (iHarding fc Mnslimov 


20nhll . 


There has also recently been a substantial increase in the populatio n of known 


VlSPs 


through discovery of new radio MSPs in Fermi unidentihed sources flAbdo et ah 


20131). 


Many of these are nearby (within 1 kpc) and a number are relatively bright, indicating that 
the existing radio surveys were incomplete (or insensitive to the detection of many MSPs). 
All of the above factors (more sources characterized by higher pair multiplicities and larger 
maximal particle energies than previously thought) make the study of MSPs as sources of 
cosmic-ray electrons and positions quite attractive. 

We have previously studied the contribution to the terrestrial electron spectrum 
by the nearby MSP PSR J0437—4715 assuming a pair-starved potential, but found the 


contribution of this nearby MSP to be negligible wit 
contribution of the much younger Geminga (see also 


lin this model. We also considered the 

hat 


Aharonian et ah 


it may contribute signi hcantly, depending on model parameters flBiisching et ah 


19951). and found 


Biisching et ah 


2 nn8ah . 


fl2008bll furthermore noted that both Geminga and PSR B0656-I-14 may 


be dominant contributors to the terrestrial positron flux, and may be responsible for an 
anisotropy of up to a few percent in this flux component. We have recently made a hrst 


attempt to carefully assess the contribution of MSPs 


clusters) to the cosmic-ray lepton spectrum at Earth fjVenter et ah 


excluding those fo und in globular 


2015|), where we have 


considered pairs originating in cascades within the magnetospheres of MSPs. However, 
since about 80% of MSPs have binary companions, in some fraction of these systems shocks 
sar winds as they interact with the companion wind or atmosphere 


may form in the pu 


(IHarding fc Gaisser 


1990 


Arons &: Tavani 


19931) . which could accelerate the pairs to higher 
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energies. It is possible that such shock acceleration occ urs in some black widow (BW) 


systems, such as PSR B1959+20 flArons fc Tavani 


19931) . Due to Fermi observations, the 


population of BWs and redbacks (RBs) has increased signihcantly. We therefore now also 
study the effect of pairs that have been reaccelerated in intrabina ry shocks of BW and 


RB systems. We furthermore i nclude Klein-Nisihna (KN) effects fjSchlickeiser fc Ruppei 


2010 


Blies fc Schlickeiser 


2 OI 2 I) when assessing the inverse Compton (IC) loss rate the 


particles suffer as they traverse the interstellar radiation held (ISRF) of the Galaxy. We 
describe the assumed source properties of MSPs by hrst discussing the central expectation 
of roughly equal numbers of electrons and positrons coming from pulsar magnetospheres 
fSection l2.ll) . after which we describe our population synthesis code used to predict the 
present-day number of MSPs as well as their location and power fSection 12.21) . We describe 
an additional BW / RB source population in Section 12.31 Moving to source spectra, we 
describe our PC pair cascade code that yields realistic pair spectra fSection l3.ip . We 
also describe the spectra injected by BW and RB systems fSection 13.21) . and motivate 
why we neglect the small contribution due to primaries (Section 13.31) . We next discuss 
our assumptions regarding the ISRF (Section 14.1|) and Galactic magnetic held strength 
fSection l4.2p . which are necessary inputs to the calculation of energy losses suhered by the 
leptons fSection 14.3p . We use this together with a prescription for particle dihusion when 
solving a transport equation (Section 14.41) to calculate the spectra at Earth (Section [5]). We 
discuss our results in Section [6], while our conclusions follow in Section [3 


2. Millisecond pulsars as sources of cosmic-ray electrons and positrons 

We hrst address the question of pair production in pulsar magnetospheres f Section [2. ip . 
specihcally as this pertains to MSPs, before describing two pulsar populations we consider 
in the rest of the paper: Galactic MSPs resulting from population synthesis modeling 


























fSection I2.2p . and BWs and RBs which may further accelerate particles flowing out of the 
MSP magnetospheres in their intrabinary shocks. ISection I2.3p . 


2.1. Pair production in pulsar magnetospheres 


Produc 


by 


Stnrrock 


ion o 


electron-positron pairs in pulsar magnetospheres, hrst proposed 


(jl97ll ). is widely considered to be critical for supplying charges to the 


magnetosphere as well as plasma for the observed coherent radio e mission. The pairs can 


be effi ciently produced in electromagnetic cascades above the PCs (jDaughertv fc Harding 
19821) by 7 rays t hat undergo conversion to electron-positron pairs by the strong magnetic 


held flErber 


19661 ). These cascades are initiated by the acceleration of primary electrons 


in strong electric helds above the neutron star surface. Curvature and IC radiation from 
these particles reaches tens of GeV, creating pairs in excited Landau states. The pairs 
lose their perpendicular momentum by emitting synchrotron radiation (SR) photons that 
create more pairs. In young pulsars with magnetic helds above 1 0^^ G, the cascades 
can produce multiplicities of 10^ — 10^ pairs per primary electron flDaughertv fc Harding 


1982 


Harding fc Mnslimov 


2nilal ). The dense pair plasma will screen the accelerating 


electric held above the gap 


llMns: 


(e.g., 


imov fc Harding 


Spitkovskv 


except in a narrow gap along the last open held lines 


20041) . Screening by pairs may provide nearly force-free conditions 


20061) throughout the magnetosphere, maintaining the narrow accelerator 


and emission gaps necessary to produce the sharp caustic y-ray peaks observed by Fermi. 
The pair plasma created by pulsars hows out of the magnetosphere along open magnetic 
held lines close to the pole and provides the radiating particles for the surrounding PWNe. 


Models of PWN 


fide Jager et af 


e require high pair multiplicit y to produce the observed SR and IG emission 


1996 


Bucciantini et af 


20111 ) 


Most MSPs, because of their very low magnetic helds, have difficulty producing 
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high-multiplicity pair cascades initiated by curvature radiation if the surface helds are 
dipolar. They are able to produce cascades from IC radia tion, but these casca des do not 


have high enough multiplicity to screen the electric helds flHarding et al 


thus assumed to have pair-starved magnetospheres flHarding et ah 


20021 ). They were 


2005h that have particle 


acceleration on all open held lines up t o high altitudes. Such magnetospheres would produce 


broad 7 -ray peaks fjVenter et al. 


200911 at earlier phase than the radio peak. However, Fermi 


detected MSPs with narrow peaks in their 7 -ray light curves arriving at later phase than 
the radio peak, very similar to those of young pulsars, implying that MSPs are somehow 


able to prod uce the high multiplie r 


held region. 


y pair cascades required to screen most of the open 


Harding fc Muslimovl (l2011alJbl ) suggested that MSPs have non-dipolar helds 


near their surface that enhance the accelerating electric helds and enable creation of more 
pairs. Introducing a generic toroidal component to the dipole held that ehect ively caused 


an ohset of the PC relative to the magnetic pole. 


Harding fc Muslimovl fj2011al) were able to 


specify the held distortion with two ohset parameters, £ and 0 o, describing the magnitude 
and azimuthal direction of the shift. Physically, £ ~ 0.1 for MSPs corresponds to the 


PC o hset caused by the sweepback near t 


held 


Deutsch 


1955 


a force-free held flSpitkovskv 


the surface. 


Dvks fc 


darding 


Harding fc Muslimov 


le light cylinder of a vacuum retarded dipole 


2004), ~ 0.2 to the PC ohset from sweepback of 


2006 ), and e > 0.2 to the PC ohset by multiple helds near 


fj2011bl ) found that for magnetic helds with £ > 0.4, 


requiring moderate surface multipole components, all known MSPs were able the produce 
pair cascades by curvature radiation. 


Aside from the requirement of held distortions to produce higher pair multiplicity 
for the 7 -ray prohles, there is evidence of a non-dipolar surface held structure in MSPs 
from the study of their X-ray emission. The thermal X-ray pulse prohles of some MSPs 
show asymmetries that require ohsets from the magnetic axis of the emitting hot spot on 
the neutron star surface in order to successfully ht the light curves. Since the emission 
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l ikely originates 


flBogdanov et al. 


rom 


2007 


C heating, it i s argued that MSPs su ch as PSR J0437—471 5 


Bogdanov 


201311 and PSR J0030+0451 flBogdanov fc Grindlav 


200911 


have either offset dipoles or offset PCs. The shift of the heated PC needed for modeling the 
light curve of PSR J0437—4715, ~ 2 km, corresponds to an offset parameter £ ~ 0.6. (In 
what follows, we will adopt values of e = 0.0, 0.2, and 0.6 in our modeling.) 

Below, we discuss two classes of MSPs which we consider to be sources of cosmic-ray 
electrons and positrons. 


2.2. Galactic synthesis model for the present-day MSP population 


We implement the results of a new study by 


Gonthier et al. 


(120151) of the population 


synthesis of radio and y-ray MSPs that lead to the present-day distribution of MSPs. This 
is assumed to be an equilibrated dist ribution within th e Galaxy whose evolution has been 


described in Section 3 of the work of 


Story et al. 


fl2007l. hereafter SGH) where the radial (p 


in cylindrical coordinates) distribution was assumed to be that of iPaczyhskil fjl990l) . with a 


radial scaling of 4.5 kpc and a scale height of 200 pc, instead of 75 pc used in t 
additi on, the supernoya kick yelocity model that was implemented was that of 


lat work. In 


Hobbs et al 


(120051) using a Maxwellian distribution with a width of 70 kms“^ (resulting in an average 
speed of 110 kms“^). The Galaxy is seeded with MSPs treated as point particles with 
ages going back to the past 12 Gyr assuming a constant birth rate of 4.5 x 10“^ MSPs per 
century as obtained in SGH. The MSPs are evolved in the Galactic potential from their 
birth location to the present time when an equilibrium distribution has been established. 


We assume that MSPs are “born” on the spin-up line with initial period Pq dependent 
on the surface magnetic field Bg, which we assume does not decay with time. We assume 
a power-law distribution for the magnetic helds. As in the case of the study of SGH, the 
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simulation prefers a power-law distribution of periods Po{Bs), with an index of a-Q, with a 
normalized distribution given by the expression 

(« + 1) Sg" 


Po (Bs) = 


3 + 1 ’ 
' min 


( 1 ) 


-Dmax -Drnii 

where Bg = i?s/(10® G) and -Bmax = 10^. We consider cib and -Bmin to be free parameters, 
which are then hxed at optimum values. In the study of SGH a preferred index of —1 was 
used. However, improved agreement with the new simulation is achieved with an index of 
ClB = —1-3. 

We assume a distribution of mass accretion birth lines, from the Edding ton critical 


mass accretion rate to about 10 ^ of the critical value, following the study by 


Lamb &: Yu 


(120051) . We parameterize the mass accretion rates with a line in the P — P diagram as was 


done in Equation (5) of SGH. The intercept of this birth line was dithered using a dithering 
parameter S. The study of SGH used a ramp distribution of 6 characterized by a linear 
function increasing with 6. We found improved agreement by uniformly dithering 6 between 
0 and 2, with the restriction that the birth period Pq > 1.3 ms. 

Recently, signihcant progress has been made in obtaini ng more realis tic pulsar 


Spitkovskvl (120061) leading to the following 


magnetosphere solutions than the retarded , vacuum d i pole (Deutsch 
electro dynamic solutions were obtained by 
prescription for the pulsar spin-down power 

^sd ~ (1 + sm a) , 


19551) . Force-free 


( 2 ) 


where p is the magnetic dipole moment, G is the rotational angular velocity, c is the speed 
of light, and a is the magnetic inclination angle relat ive to th e pulsa r’s rotational axis. 


Gonsidering accelerating helds and force-free solutions. 


Lietaf 


( 120121 ) constructed solutions 


of magnetospheres Idled with resistive plasma, arriving at a very similar spin-down formula. 


Gontopoulos et al. 


(120141) considered the ideal force-free magnetosphere everywhere except 
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within an equatorial current layer, and also arrived at a similar prescription for L^d- These 
results encourage us to implement such a spin-down model into our population synthesis 
code. Using a dipole moment of / 2, where R is the stellar radius and Bg the 

surface field at the pole, and equating Lgd to the rate of rotational energy loss yields the 
expression 

^2 ^ Gc^IPP 
® i?® (l -|- sin^ a) 

Integrating this equation over the age t of the pulsar provides the expression for obtaining 
the present-day period P(t) 


/’^ = ^o^ + ^(l + sin^«) B^t. 


(4) 


We assume i? = 12 km and MSP mass Mmsp = 1-6 Mr?, , where is the mass of 


the Sun. We use the prescription outlined in Section 2 of 


Pierbattista et al 


(1201211 to 


obtain the moment of inertia, which with these values of R and Mmsp yields a value of 
J = 1.7 X 10^^ gcm^. While there is growing evidence that the inclination angle becomes 


aligned with the neu 


(jjohnston et al. 


20071: 


ron star’s rotation al axis with time in the case of normal pulsars 


Young et al.l 


2010 ), we do not consider such an alignment model in 


the case of MSPs. 


Figure [T] indicates histograms of period logio(-P )5 period derivative logio(P’), surface 
magnetic field log^o(-^s), and distance d characterizing the simulated present-day Galactic 
MSP population. Figure |2] shows several best-£t simulated and observed radio properties 
(log;^o(-P)) logio(-P)) characteristic age logio(Tc), and logio(i?s)) of radio-loud MSPs detected 
in 12 radio surveys. The output from this simulation predicts the location as well as P and 
P of roughly 50,000 Galactic MSPs, which we use as discrete sources of relativistic electrons 
and positrons in the calculations that follow. 
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2.3. MSPs in binary systems — BWs and RBs 

The majority of MSPs (about 80%) are in binary systems, and a subset of these, the 
BWs and RBs, may contain strong intrabinary shocks that can further accelerate the pairs. 
BWs are close binary systems, with orbital periods of hours, containing a rotation-powered 
MSP and a compact companion having very low mass, ~ 0.01 — 0.05 Mq. The companion 
stars in BWs undergo intense heating of their atmospheres by the MSP wind, which drives 
a stellar wind and rapid mass loss from the star. A shock will form in the pulsar wind at the 


IHardine: & Gaisser 

1990; 

Arons & Tavani 

1993) 


19931) . RBs are si milar systems 


, except that the 


2011 ). The MSPs in 


companions have somewhat higher masses, ~ 0.1 — 0.4 Mq fiRobertsI 
both types of system are typically energetic, with Lgd ~ 10^^ — lO^^ergs”^. Figure [3]is a 
schematic view of a shock formed between the colliding pulsar and companion star winds. 

Before the launch of Fermi these systems were rare, with only three BWs and one 
RB known. The large amount of material blown off from the companion stars absorbs 
and scatters the radio pulsations from the MSPs, making them difficult to detect at 
radio wa velengths. In th e last few years, radio searches of Fermi unidentihed y-ray point 


sources (iRav et ah 


20121) have discovered 14 new BWs and 6 new RBs to date, making a 


present total of 24 of these systems. In order to assess the contribution of these systems 
to the Galactic cosmic-ray positrons, we compiled a list of public detections, plus some 
measured and derived quantities (see Tables [T] and [2]) . In deriving the spin-down luminosity 
and surface magnetic helds for the pulsars in these systems, we used an MSP radius of 
R = 9.9 X 10® cm and moment of inertia of 1.56 x 10"^® gcm^, in order to be consistent with 
our pair cascade model assumptions (Section l3.ip . 

Evolution models and pop ulation synthesis of MSP binary systems yield a birthrate for 
BW systems ~ 1.3 x 10“^yr“^ flKiel fc Taamll2013l ). Taking an age of the Galaxy around 12 
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billion years, there may be a total population of several thousand BW systems. Since only a 
small fraction of these have been discovered, it is harder to estimate how many undiscovered 
BW and RB systems are within several kpc of Earth. Conservatively, the known nearby 
population may be ~ 10% of the total, or around several hundred. By considering only the 
24 known BWs and RBs, we are obtaining a lower limit to the cosmic-ray flux contribution 
by binary MSPs. 


3. Models for pair injection spectra 

3.1. Compntation of pair spectra from pulsar polar caps 

We calculate the spectra of pairs leaving the MSP magnetosphere using a code that 
follows the development of a PC electron- positron pair cascade in the p ulsar magnetosphere 


(details of the calculation can be found in 


Harding &: Muslimov 


2011b|) . The pair cascade is 


initiated by curvature radiation of electrons accelerated above the PCs by a parallel electric 


held, derived assuming space-charge-limitec 

how (i.e., free emission of particles from the 

neutron star surface; 

Arons & Scharlema,nn 

1971 

1). A fraction of the curvature photons 

undergo magnetic pair attenuation ( 

ErbeJ 

1966; 

Daughertv & Harding 

1983 

), producing 


a first-generation pair spectrum which then radiates SR photons that produce further 
generations of pairs. The total cascade multiplicity M_|_ (average number of pairs spawned 
by each primary lepton) is a strong function of pulsar period P and surface magnetic held 
strength Bg, so that many pulsars with low magnetic helds and long periods produce either 
few or no pairs for dipole held structure (e: = 0), leading to a pair death line in the PP 
diagram. 


However, as discussed in Section 12.11 the sweepback of magnetic held lines near the 
light cylinder (where the corotation speed equals the speed of light) as well as asymmetric 
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currents within the neutron star may cause the magnetic PCs to be offset from the dipole 


axis. W e adopt the distorted magnetic field structure introduced by 


Harding fc Muslimov 


f)2011bl) that leads to enha nced local electric field s 


below the pair death line. 


, boost ing pair formation, even for pulsars 


Harding fc Muslimov! fj2011bll considered two configurations for 


the dipole offset in which the magnetic field is either symmetric or asymmetric with respect 
to the dipole axis. Sweepback of the global field would pro duce asymmetri c offsets, while 


the observed offset in the MSP J0437—4715 is symmetric flBogdanov 


2013 1. We adopt a 


symmetric field structure for calculating the pair spectra of MSPs in this paper. In the 
symmetric case, the magnetic field in spherical polar coordinates {rj, 6, 0) is 


B 




T COS 9 + - 6 {1 + a) sin 9 — <p e sin 9 cos 9 sin((;/) — ^q) 


(5) 


where Bg is the surface magnetic field strength at the magnetic pole, rj = r/R is the 
dimensionless radial coordinate in units of neutron star radius R, a = e cos(0 — 0o) is 
the parameter characterizing the distortion of polar field lines, and 0o is the magnetic 
azimuthal angle defin i ng the meridional plane of the offset PC. Using this field structure. 


Harding &: Muslimovi (1201 Ibl) derive the component of the electric field parallel to the local 


magnetic field, E», t 


r at acce lerates electrons. We have used the E» of Equation (11) of 


Harding &: Muslimovi fj2011b|l that corresponds to a symmetric offset and use these field 


structures to accelerate the electrons above the PC to simulate the pair cascades. The pair 
spectra (Figure 0]) are characterized by P, P (or equivalently, Bg via Equation |3]), and 
offset parameter £. From our simulations, we find that about ~ 1% of Lgd is tapped to 
generate the pairs. 


We used a grid in P and Bg encompassing P = (1,1.8, 2, 2.5, 3,4, 5, 7,10, 20, 50,100) ms, 
and Bs = (1,1.5, 2, 3, 5, 8,10,15, 20, 50). For each source in the present-day MSP population 
with predicted values of P and P fSection l2.2p . we found its associated pair spectrum 
by interpolating spectra on this grid. We used an inclination angle of a = 45°, mass 
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^MSP = 2.15Mq, radius R = 9.9 km, and moment of inertia I = 1.56 x lO"^® gcm^ for all 
MSPs. We adopted an equation of state with larger Mmsp here (and associated smaller 


I) compared to that used in the 


measured masses Mmsp ~ ^Mq (IDemorest et ah 


population coc 


e (Se ction 12 .2p . since some MSPs have 


2010 l) . and this enhances pair multiplicity. 


However, this discrepancy is removed by considering a large range of £, since the latter 
simulates a large range of pair multiplicities that would correspond to different equations of 
state, and thus different values of Mmsp- We used dipole offsets of £ = (0.0, 0.2, 0.6) and set 
00 = 7r/2 (this parameter controls the direction of offset of the PC). 


We use the above spectra as input for the calculation of the positron component from 
the population-synthesis sources fSections 12.21 and 151). Since MSPs are not surrounded by 
nebulae that can trap the pairs and degrade their energy before escape, we can assume that 
the pair spectra emerging from the MSPs are good representations of the actual source 
spectra. 


3.2. Spectra from particles accelerated in the intrabinary shocks of BWs and 

RBs 


We assume that the pairs escaping from the pulsar magnetosphere may be further 
accelerated in the intrabinary shock that originates between the pulsar and companion 
winds in BW and RB systems. Acceleration of leptons at a large distance outside the 
pulsar light cylinder is necessary to account for the extended SR emission observed from 
PWNe. Such acceleration is thought to occur at or near the termination shock in the 


pulsar wind (iKennel fc Coroniti 


1984J ) that is conhned by the sub-relativistic expansion 


of the surrounding supernova shell. The acceleration mechanism near the pulsar wind 
termination shocks is not understood, but is known to be highly efficient, since the 
bolometric luminosity of the Crab nebula is about 20% of the pulsar spin-down luminosity 
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and the inferred maximum pa rticle energy, 


voltage across open field lines fide Jager et al. 


10^^ eV, is at least 10% of the available 


19961) . The pulsar wind termination shock is 


relativistic and perpendicular, so that the diffusive hrst-order Fermi mechanism becomes 
problematic u nless most of the magnetic energy is converted into particle energy upstream 


fiSironi fc Spitkovskv 


of the shock fiSironi fc Snitkovskv 

3 


2nilal ). However, either s hock-driven reconnect ion 


2011bf) or strong electromagnetic waves (lAmano fc Kirk 


3 ) 


20131) could 


cause demagnetization, enabling diffusive acceleration to proceed. 

Regardless of the acceleration mechan ism, the maxim um particle energy will be limited 


by the universal scaling, 


vBRs/c flHardine 


19901), where n is a bulk flow velocity. 


B is the magnetic held strength, and Rs is a scale size of the system. In the case of shock 
acceleration, the maximum energy comes from a balance between the minimum acceleration 
timescale, set by the particle diffusion, and the timescale for escape from the shock of radius 
Rg. However, for leptons, the timescale for SR losses is shorter than the escape time and 
the maximum energy will be set by balancing the acceleration timescale with the SR loss 
timescale. 

We assume that the reaccelerated shock-accelerated spectrum will be an exponentially 
cut off power law with spectral index of —2 

Eq 


Qi{E) = Qo^iEq exp 


E, 


cut 


( 6 ) 


with the index i indicating the source, Qo,i the normalization factor, and Eq the particle 
energy at the source position. In order to estimate the maximum (cutoff) energy, we 
balance the energy gain rate from shock acceleration, assuming Bohm diffusion, and the SR 
loss rate that particles experienc e in the strong magnetic held at the shock radius. This 


199nh 


leads to the following expression flHarding fc Gaisserl 

■ 3 (^- 1 )^'/' 


E, 


cut 




U(e + i)J 


TeV, 


(7) 
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with Pjns the pulsar period in milliseconds, an = a/(10^^ cm) the binary separation, 
and £ the shock compres sion ratio. This is slightly different from Equation (34) in 


Harding &: Gaisserl fjl990l) . since they assumed that the shock distance from the pulsar is 


Ts ~ a — i?*, with i?* the companion radius. For BWs and RBs, the shock is close to the 
companion star, and we assume Vg ^ a, leading to the modihed expression given above. The 
binary separation may be found as follows (given the extremely small eccentricities of these 
systems) 


a = 


G (Mmsp + M 


comp ) 




1/3 


Pv 


2/3 


b > 


( 8 ) 


with Mmsp the MSP mass, Mcomp the companion mass, Pb the binary period, and G the 
gravitational constant. We have listed the inferred values of an and Pcut for each of the 
detected BWs and RBs in Tables [T] and |2j 


Now, we can normalize the spectrum (e.g.. 


Biisching et al. 


2008b|) using 


Qi (IEq = [M+(P, Ps, e) + 1] nGj(-P, 


QiE{) dE() '^p,max-^sd; 


(9) 

( 10 ) 


with M_|_ the pair multiplicity, r/p,max the efficiency of conversion of spin-down power 
Lsd = dvr^/PP”^ to particle power (or shock efficiency), and hcji^P, Bs,e) the Goldreich- 
Julian particle outflow rate , appropriate for offset-dipole helds (see Equation [3] of 


Harding &: Muslimov 


201 Ibh characterized by an offset pararneter g fSection 13.ip . The 


latter is similar to the classical expression (iGoldreich fc Julian 


196911 


’^GJ = 


2cHpcPGJ dvr^PsP^ 


e 2ceP^ ’ 

with A-pc the area of one PG, and pcj the Goldreich-Julian charge density. In Equation (C 
we therefore normalize the spectrum to the total (primary plus secondary) current. We 


( 11 ) 


found bv interpolating values on a grid o 

' P and P„, while we calculated iiaAe) directlv 

from the cascade code ( 

Harding & Muslimov 

2nilb 

)• 
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The above is a system of two equations and two unknowns, Qo,i and -Emin (when hxing 
hp.max)- We hnd that the spectrum of Equation (jH]) can only be normalized for some choices 
of M+, Ecut, and ?7p,max- Figure 0 shows contour plots of log^g (-Emm/Ecut) vs. logio(M-^) 
and logio(-^cut) assuming Pms = 3, Bg = 5, R = 9.9 x 10^ cm, and / = 1.56 x 10^^ gcm^. 
Panel (a) is for r/p max = 0.1, while panel (b) is for r/p^max = 0.3. Values near unity (dark red 
regions, i.e., the lower left corners) indicate that no solution could be found for the given 
parameters. Fixing r/p^max, one can see that for a hxed value of -Ecut, some minimum value 
of M+ is required in order to hnd a physical solution -Emin < Ecut- This is because a higher 
M+ will raise Qo^i, allowing Equation (E]) to be satished. A higher value of -Emin has the 
same effect. For an even higher value of M+ (typically associated with a higher value for 
e) than the critical one needed to hnd a physical solution, the constraint on -Emin relaxes, 
and one hnds a smaller ratio -Emm/Ecut, and therefore a spectrum spanning a larger energy 
range. In other words, if is too low for a hxed value of -Ecut, h is not possible to satisfy 
the constraint of the total power (Equation [lO]). To solve this problem, we decreased ?7p,max 
systematically until we found a solution. Comparison of panel (a) and panel (b) indicates 
that a smaller value of T^p^max will relax the power constraint, so that solutions may be found 
for larger regions in (M+,.Ecut) space. 

With the solutions of source spectra in hand for the population of 24 BWs and RBs 
considered (Tables [1] and |2]), we may next calculate their transport through the Galaxy 
(Section 03]). 


3.3. Neglecting the primary component from population-synthesis MSPs 

We have noted that the secondary component almost always vastly dominates the 
primary component in the case of the BWs / RBs fSection 13.21) . i.e., usually M+ S> 1. This 
is due to the fact that multiplicities grow very rapidly with e. Even in the case of £ = 0.0, 
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while the primary spectra may dominate the secondary spectra for some low--Bs and large-P 
pnlsars (which wonld imply M_|_ <C 1 ), there will always be pnlsars with high enongh 
and short P (i.e., M_|_ ~ 100 — 1000) so that their secondary spectra will dominate the 
cnmulative flnx contribntion from a popnlation of pnlsars. This means that the cnmulative 
spectrnm from the BW and RB pnlsars will be dominated by secondary, and not by primary 
spectra^ 

On the other hand, for the MSPs from onr popnlation synthesis model, where we 
assnme no shock acceleration, the primaries may form nearly mono-energetic spectra at 
very high Lorentz factors 7 ~ 10^“®, depending on held-line cnrvatnre, i.e., colatitnde, and 
also P and Given this small energy range (the spectrnm is almost a 5-distribntion), one 
might think that this component may leave a distinct signatnre in the total spectrnm of 
particles leaving the pnlsar magnetosphere. However, when combining primary spectra from 
several pnlsars, and following their transport throngh the Galaxy to Earth, the cnmnlative 
primary spectrnm will have been smeared ont due to the different source locations and 
properties. The primary spectra should also be at a lower intensity than the secondaries, 
given the typical multiplicities encountered for the Bg and P values of the closest MSPs. 
Furthermore, if there would have been any signature at high energies ~ 10 TeV, where the 
secondary spectra drop off in this case, this will be completely masked by the contribution 
of the BW / RB. 

Given the above arguments, we do not include the primary spectra from the MSP 
synthesis population since they should not have an impact on our results. 

^Neglecting the primary spectra in this case would imply setting -|- 1 M_|_ when 

solving Equation (I9l). While we have not done this, the effect would be negligible, given the 
large values of M_|_ in some cases, even for e: = 0 . 0 . 
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4. Galactic transport of injected leptons 


4.1. Interstellar Radiation Field 


Knowledge of the spectral and spatial properties of Galactic ‘background photons’ is 
important for calculations of IC losses suffered by leptons propagating through our Galaxy. 
The relevant photons are optical ones produced by the population of stars in the Galaxy, 
in addition to infrared (IR) photons that are the result of sc attering, absorption , and 


re-emission of the stellar photons by dust in the ISM; see e.g., 


Porter et ah 


(1200811 ■ 


The GALPROP code (IStrong fc Moskalenko 


1998f) includes a detailed model for 


this ISRF that incorporates a stellar population model (i.e., a luminosity distribution 
derived from 87 stellar / spectral classes distributed in 7 geometric locations within the 
Galaxy), dust grain abundance and size distribution models, as well as the absorption 
and scattering efficiencies of the latter which enable radiative transport calculations for 
stellar photons propagating through the ISM. While the ISRF is inherently anisotropic and 
inhomogeneous, with the bulk of the photons leaving the inner Galaxy, the ISRF model 
used by G ALPROP assumes azimu t hal symmetry and a cylin drical geometry. For more 


details, see 


Moskalenko et ah 


fennbh : 


Porter et al. 


(12996 


29981) and references therein. 


For our purposes, we only need average photon energy densities to calculate the 
total IG loss rateci, since this is the quantity needed to solve the transport equation 
(see Equation |22]). We find that the GALPROP ISRF is adequately approximated by 


three black 
Figure 2 of 


3odv components 


Venter et al. 


(optical, IR, and cosmic microwave bac kground or GMB; see 


2915h . We follow iBlies fc Schlickeiserl (129121) in distinguishing two 


main spatial regions: the Galactic Disk and the Galactic Halo. For the Disk, we use their 


^However, the temperature Tj of each blackbody component j is needed when implement¬ 
ing KN corrections; see Equation (fT6|l . 
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values of Un ^t = Um = 0-4 eVcm a nd Uqub = 0.23 eVcm which is similar to the 


values of ISchlickeiser Sz Ruppell fj2010j) . while for the Halo, we use f/opt = 0.8 eVcm 
f/iR = 0.05 eVcm ~^, and Ucmb = 0 . 23 eV cm“^. The dust is assumed to follow the Galactic 


gas distribution (iMoskalenko et ah 


20061 ) ■ which tapers off strongly with perpendicular 


distance above the Galactic Plane, leading to less absorption of optical photons (and 
therefore a larger value for Uopt and a reduced value of Gir in the Halo). 


4.2. The Galactic Magnetic Field 


Had (120091) noted that there are hve observational tracers of the Galactic magnetic 
held. These are polarization of starlight (indicating that the local held is parallel to the 
Galactic Plane and follows the local spiral arms); polarized thermal dust emission from 
molecular clouds (indicating held enhancement upon cloud formation via compression of 
the ISM, and that magnetic helds in these clouds seem to be preferentially parallel to the 
Galactic Plane); Zeeman splitting of spectral emission or absorption lines from molecular 
clouds or from OH masers associated with HII or star forming regions (indicating large-scale 
reversals in the sign of the line-of-sight component of the median held, and that interstellar 
magnetic helds are apparently preserved through the cloud and star formation processes); 
dihuse SR radio emission (used to estimate the total and ordered or regular held strength); 
and Faraday rotation of linearly polarized radiation from pulsars and extragalactic radio 
sources (giving a measure of strength and orientation of the line-of-sight component of the 
magnetic held). The combination of the latter with measurements of total intensity and the 
polarization vectors (from SR) all ows one to distinguish between three held components: 


regular, anisotropic, and random (IBeck 


20091). 


For our purposes, we are interested in an average held strength that would determine 
the SR loss rate (Equation [13] below), and not so much in the overall Galactic held 
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structural (which is still under debate). The total held has been estimated to be around 
6 /xG, averaged over a distance of 1 kpc around the Sun, using SR measurements and 
equipartition arguments where the magnetic energy density is set equal to that of cos mic 
rays. This number incr eases to 10 pG closer to the inner Galaxy (see 


and references therein). 


Han et al. 


Beck 


20091. 


fj2006h used a combination of dispersion and rotation 


measures of over 500 pulsars and found that the regular magnetic held decreases from 


~ 6 fiG near a Galactocentric c 
near the Sun (see Figure 11 of 


istance of 2 kpc to ~ 1 /iG near 9 kpc; the value is ~ 2 /iG 


Han et al. 


20061) . The latter should be compared to recent 


measurements o 


(iBurlaga fc Nest 


the i nterstellar magnetic held by Voyager 2 which yielded 3.8 — 5.9 fiG 


20141) . Furthermor e, the mean regu lar held as function of latitude is 


inferred to vary between ~ ±5 fiG flHan et al. 


20061) . Fields in interarm regions are 


seemingly weaker than those in spiral arms . For example, 


Norma arm was found to be 4.4 ± 0.9 /iG flHan et al 


lie average regular held in the 


20021) . The regular magnetic held 


has only a weak vertical co mponent of = 0 . 2 — 0.3 mG, directed from the southern to 


the northern Galactic Pole (iHan fc Qiao 


19941 ). lOrlando fc Strong (120131 ) inferred values of 


~ 2 fiG, ~ 5 /iG, and ~ 2 fiG for the local regular, random, and anisotropic held components 
in the Disk via Galactic SR modeling. The average total held, ho wever, decrease s when 


taking into account its rapid decay with height above the Plane. 


Delahave et al. 


(12010 


hereafter DIO) argue that SR losses depend on the mean of the squared magnetic held, so 
that one should include all components in the following way: 


bsr = + <-®2> + (A">. 


( 12 ) 


with Hr the regular held, the irregular, anisotropic held aligned with the regular on e, and 


Hi the isotropic or random held; also, (R^) = (Rr)^- Results from 


Jaffe et al. 


(l2010[) lead 


^Kistler et al. 


(120121) raised the additional issue of particle transport in a turbulent mag¬ 


netic held, which we will briehy consider in Section [61 
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to values of up to i?sR ~ 6 /iG for fields in the Galactic Disk. However, if an exponential 
decay function for the vertical behaviour of the magnetic held is assumed, and Hsr is 
averaged over a spherical volume of radius 2 kpc, DIO hnally obtains an average value of 
Bsr ~ 1 — 3 /iG. 


4.3. Total Leptonic Energy Loss Rate 


The SR loss rate is given by 


Esr — 


AaT'cUBE'^ 


^2 ’ 


3 {rrieC^) 

with rUe the electron mass, E the particle energy, and Ub the magnetic energy density 

j^2 

Ub = -rr— = 0.0986n eV cm“^ 


(13) 


Stt 


(14) 


for a Galactic held of 62 = B/{2.0 /iG). The general expression (including KN ehects) for the 
IG loss rate (for target photons of energy density Uj, and j signifying diherent blackbody 
components associated with temperature s Tj) may be approximated as (for details, see the 


Appendix of 


Schlickeiser fc Ruppelll2010l) 


^icj - 


AarcUiE^ 


7kn,- 


3(mec2)^ 7 knj + 7^’ 
with 7 the particle Lorentz factor, and the critical KN Lorentz factor dehned as 

S'/EirieC^ 0.27meC^ 

7knj =-~- 


(15) 


(16) 


Stt kBTj kBTj 

The IG loss rate for particles with Lorentz factors above 7 kn,j is seve rely suppressed. If 7 
x 2 ' we recover the well-known expression for the Thomson limit flBlumenthal fc Gould 


1970) 


• _ AaTcUjE 

-^ICJ — 


3 {meC^) 


2 • 


(17) 
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By considering varions Galactic soft-photon target fields (IR, CMB, and optical) with 
respective energy densities Uj, we note that the KN correction is only necessary for optical 
photons, where 7KN,opt ~ 10^. 


Previously (jVenter et ah 


2 OI 5 I) . we assumed that we could approximate all losses as 


being in the Thomson limit for all cases. Since all loss terms (SR and IC, for the different 
soft-photon components) have the same functional dependence on energy in this case, 

E oc E'^U, where U can indicate either magnetic or soft-photon energy density, we could 
define one single loss term using an effective magnetic field ReS that takes into account 
both SR and IC losses. We found a value of ReS ~ 7 /rG for both the Plane and the Halo, 
given the typical values used for Uj and R. These Reg values are the same in both regions 
because Ropt goes from 0.4 eV cm“^ in the Plane to 0.8 eV cm“^ in the Halo, while Um goes 
from 0.4 eVcm“^ in the Plane to 0.05 eV cm~^ in the Halo. In addi tion , the actual magneti c 


held drops from R ~ 3 /rG to R ~ 1 /rG f Blies fc Schlickeiser 


2 OI 2 I). In 


Venter et af 


( 1201 ^, 


however, we decided to use a slightly lower value of Reg = 5 /rG in view of the fact that the 
Thomson limit would overestimate the losses. 

For this paper, we introduce two loss terms, thereby separating those in the Thomson 
limit, and the one in the KN limit. We denote this as follows: 


T^total — I E' 


^Thom 


+ -EkN — ( -^SR + E] 


4c, IR 


+ E 


IC,CMB 


+ E 


IC,opt • 


(18) 


To calculate RThom, we formally set Ropt = 0 eVcm We can combine the rest of the 
terms into one by dehning an effective magnetic held, Reg, as before: 

4(JtcR^ 


-E^Thom — Esb. + RiC.IR + 77lc,CMB — —- 

3 


\Ub + Rir + Rcmb] 


R = 


3 {nieC^y 


4c 


9 V rUeC^ 


= 1.58 X 10 


-15 I EeS 


1 fiG 


(19) 

( 20 ) 
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We find values of = 5.2 — 5.9/iG in the Plane, and i?efr = 3.6 — 4.6/iG in the Halo 
(the range stemming from the fact that we considered B to be in the range 1 — 3/iG; see 

Section 02]) • 


We have to treat the optical component separate l y, sinc e the KN effect will become 


important in this case. Following 


Blies fc Schlickei.serl (l2012h . we consider only the most 


dominant optical component, which we approximate by a black body with Topt = 5, 000 K, 
replacing Uj by f/opt and 7KN,i by yxN.opt in Equation ([15]): 

'' TTCGoptE^ 7KN,opt 


Ekn — 


( 21 ) 


3 (mec2) 7KN,opt + 7^' 

The particles will traverse r egions having different ( line-of-sight-averaged) f/opt. Furthermore, 


the optical stellar model of 


Wainscoat et af 


(119921) used to calculate the ISRF gives a scale 


height of 0.27 — 0.325 kpc for the stars of Topt ~ 5, 000 K, while the scale heig 


it for the 


Galactic magnetic field varies between ~ 0.1 — 4 kpc (DIO; 


Orlando fc Strong 


201311 . In 


view of the uncertainties associated with obtaining a line-of-sight-averaged B^q and f/opt 
for each source, and given the fact that our transport model considers only one spatial 
dimension, in what follows we consider two extreme cases of minimal and maximal total 
losses Ttotai to bracket our particle flux results: (1) Tefr = 3.6 /iG and f/opt = 0.4 eV cm“^; 


and (2) Teff = 5.9 /iG and f/opt = 0.8 eV 


cm 


-3 


4.4. Solution of the transport eqnation 


In order to transport the pairs created in the MSP magnetospheres to Earth, we use 
the following Fokker-Planck-type equation that includes spatial diffusion and energy losses: 

f)r) f) / ' \ 

^ = V • (/C ■ Vue) - — (Ttotaine) + (22) 


with rig the lepton density (per energy interval). Also, /C denotes the diffusion tensor and 
-Etotai the total energy losses, while S is the source term. 
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Since MSPs are quite old (ages of ~ 10^*^ yr), and have very small time derivatives of 
their period P, we assume a steady-state scenario {d/dt = 0). We furthermore assume a 
uniform ISM, and thus invoke spherical symmetry such that Ug only depends on distance 
r = d from Earth. For a scalar diffusion coefficient k, Equation fl22|) now reduces to 


0 = 1A (r\^ 

r 2 Qr V dr 


d 

( -^total^e ] Q' 


(23) 


We incorporate the energy losses as explained in Section 14.31 and assume that the 
diffusion coefficient is spatially independent (so that /C becomes a function of energy only, 
k{E)) and we as sume a power law energy dependence and (as motivated by quasi-linear 


theory; see, e.g.. 


Maurin et ah 


200211 


k{E) = Kq 


QD 


(24) 


We assume tvpical values of an = 0.6, = 

1 GeV, and Kn 

= 0. 

kpc^Mvr ^ ~ 

3 X 10^® cm^s ^ fe.g.. 

Moskalenko & Strong 

1998; 

Malvshev et ah 

2009 

Grasso et ah 

2009 


Feng fc Zhang 


20151) . The value for Kn is an indication of the efficiency of the diffusion 


process at a particular energy (e.g.. 


Maurin et af 


20021 ). while ctd is inferred from the 


measured ratio o 
the Galaxy (e.g.. 


boron to c arbon and character izes the escape time of cosmic rays from 


Blasi 


2009 


Genolini et af 


201511 


For the source term, we consider ~ 5 x 10^ Galactic MSPs from the population 
synthesis code (Section 122]), and = 24 for the BW / RB case fSection l273|l . For the 
pulsar in our synthesis population, we assign a pair spectrum Qi{P, Bs,e, E), as calculated 
in Section Ism for the corresponding simulated values of P, Bg, and £. We model this as 

N 


S = ^Q.(P,B.,B)i(r-r„,i). 


(25) 


Here, tq,, are the source positions. For an infinite system. Equation (l22ll is solved by the 
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following Green’s function (e.g., DIO; 


Blies&Schlickeiser 


201211 : 


G(r, ro, E, Eq) = 


|r - rol 

A 


-E'total (^A) 

with Eq the particle energy at the source, and the square of the propagation scale is 
characterized by 

rEo 


(26) 


X{E, Eq 


= Ao 


E Etotal{E') 


dE', 


An — 


1 

Eq V E, 
4ki 


En 


au 


E \E, 


E 


od 


0 


(27) 

(28) 
(29) 


(cxd — 1) &0 ’ 

and Q{Eq — E) the Heaviside function. The latter is used to ensure that A > 0. The lepton 
flux may then be found using 

c 


0e(r,^) = 


dvr 



G(r, ro, E, Eq)S dEod'^VQ. 


(30) 


While the finite boundary of the Galactic Halo should impact the solution, this effect is not 
too large for GeV leptons, for which the propagation scale is only a few kpc (DIO), and we 
neglect it here for simplicity. Our results will indicate that our predicted MSP contribution 


becomes si gnificant above ~ 1C 


neglected fjStrauss fc Potgieter 


GeV , so that 


201411 . Indeed, 


he effect of so l ar mo dulation may safely be 


Accardo et ah 


(2014) noted t 


has no effect on the newly measured PF by AMS—02, although 


lat modulation 


Aguilar et ah 


f 2014l l claimed 


that they see the effects of solar modulation up to ~ 10 GeV in their electron and positron 
data. 


We found that in order to have smooth output spectra, we had to treat nearby and 
distant sources separately. We used a logarithmic grid for the particle source energies 
Eq of the distant sources {d > 1 kpc), which is strongly refined in Eq for the nearby 
(d < 1 kpc) ones (~ 100 sources) when solving Equation fl30ll . This was necessary, since 
there are “poles” in the Green’s function when E Eq, and A ~ 0, so G(r, fq, E, Eq) —)■ oo 
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(Equation [26]). These singularities are however, removable, in the sense that a very hne 
grid in Eq results in a hnite integrand for Equation flSU]) . while the Heaviside function 
formally avoids Eq = E. 

Figure |6] shows a comparison of transported spectra involving the Galactic synthesis 
MSP component, in the Thomson limit (including all background photons, plus SR losses) 
vs. the KN limit (i.e., SR, Thomson limit for IR and CMB, but KN limit for the optical 
photons). See Section 14.31 for details. We compare Disk and Halo scenarios. In the Thomson 
limit, these imply the same value of ReS = 7pG (indicated by solid lines; different values 
for £ are distinguished by the colors); however, in the KN limit, we have to separate the 
optical photon component, and we indicate the relevant values in the legend (dashed lines 
are for the Halo, for Refr = 3.6 pG and Gopt = 0.8 eVcm“^, while dotted lines indicate 
Reff = 5.9/iG and Ropt = 0.4 eVcm“^, for the Disk). It is noticable that the particles at 
higher energies suffer fewer losses in the KN case, given the reduction of the losses above 
~ 160 GeV in this regime, and hence this raises the transported spectrum somewhat. The 
largest enhancement of particle flux occurs for the Halo case, given the low value of Reff. 
We also show the effect of changing the normalization of the diffusion coefficient. For a 
larger kq (cool colors), the flux is lower, while the opposite occurs for a smaller value of 
Ko (warm colors). One may view the latter case as a pile-up of particles, and one can also 
observe a transfer of high-energy particles to lower energies, given the slower diffusion, as 
evidenced by the change in slope at lower energies. Lastly, an increase in flux with e is 
evident, given the larger value of M+ implied by a larger value of £. 


5. Results 

Figure [7] shows the synthesis and BW / RB spectra transported to Earth, as well as 
the sum of the synthesis and BW / RB components (see legend). Dashed lines are for 


30 


Beff = 3.6/iG and f/opt = 0.4 eVcm“^, while dotted lines indicate i?eff = 5.9/iG and f/opt 
= 0.8 eVcm”^ (i.e., minimal and maximal particle losses). The different valnes for e are 
indicated by different colors as noted in the legend, and we assnmed Kq = 0.1 kpc^ Myr“^ 
and ?7p,max = 0.1. Figure [8] is the same, but for r^p^max = 0.3. While the synthesis component 
contributes mostly at tens of GeV, the BW / RB component contributes at thousands of 
GeV. It is also noticable that the BW / RB contribution is higher in this case due to the 
larger maximum shock efficiency. The discontinuous jump in the total spectrum is caused 
by the fact that the BW / RB component is the sum of a small number of spectra that 
cut off at particular values of the minimum particle energy i^min (calculated in each case by 
suitable normalization of the various binary injection spectra; see Section [3.21) . and that 
the flux of this component dominates over that of the synthesis component (which cuts 
off around ~ 30 GeV), making these low-energy cutoffs more evident. Furthermore, the 
spectral variations at low energies for the BW / RB component may be attributed to the 
fact that we are adding only 24 detected sources to obtain (a lower limit of) the cumulative 
contribution of binary MSPs to the cosmic-ray flux. Such variations should be smoothed 
out if a larger number of sources is used in this calculation. 

We next investigated the effect of varying the energy dependence of the diffusion 
coefficient by varying the parameter od (see Equation [21] and Figure [HD. We fixed 
Reff = 3.6/iG, f/opt = 0.4 eVcm“^, also setting r^p^max = 0.1, Kq = 0.1 kpc^Myr”^, and 
choosing values of ctD = 0.3 and 0.6, given the uncertainty of this parameter. (In the rest 
of the paper, we fix od = 0.6, unless stated otherwise.) Different values of £ are indicated 
in the legend, as in previous plots. A smaller value of ctD corresponds to relatively lower 
diffusion coefficients above the break energy of Rnorm =1 GeV. This has the same effect as 
assuming a smaller normalization kq, i.e., the spectra at Earth are relatively higher due to 
increased particle density. This effect is even more evident at higher energies, where the 
BW / RB component dominates, leading to significant uncertainties in this component’s 
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flux. The situation is opposite for source particle energies smaller than T^norm, and one can 
see the opposite effect at terrestrial particle energies below ~ 200 MeV. 


Figure 


indicates the “backgr ound” electron and positron fluxes predicted by 


UtAnmur iiviac 

Ackermann et ah 
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2012 
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Adriani et ah 
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Aguilar et ah 


2014J) . We indicate synthesis spectra plus BW / RB spectra (we assume equal numbers 


of positrons and electrons) for dipole offsets of £ = (0.0, 0.2, 0.6) and combinations of 
(Reff, Ropt) = (3.6/iG, 0.4 eVcm“^) and (ReS, Gopt) = (5.9/xG, 0.8 eVcm“^), i.e., minimal 
and maximal losses. The various curves are distinguished in the Figure caption. We set 
Kq = 0.1 kpc^Myr"^ r]p^uiax = 0.1. 


Figure [TTl indicates the case for kq = 0.1 kpc^ Myr“^, = 0.3, while Figure O and 

Figure fT3] are for Kq = 0.01 kpc^ Myr“^, r/p^max = 0.1 and Kq = 0.01 kpc^ Myr“^, r^p^max = 0.3, 
respectively. As before, we see that the BW / RB contribution is higher for a larger shock 
efficiency, and that all components are higher for a smaller diffusion coefficient. This is due 
to a pile-up effect which boosts the particle density. We indicate the effect of changing ajy 
in Figure [HI where one can see that the flux increases for smaller values of od, as noted 
earlier. For comparison, we show in Figure [15] the results when using the background model 
of DIO (we use their secondary positron flux as well as the sum of their secondary electron 
flux and primary electron flux originating in distant SNRs, as indicated in their Figure 14), 
where we assume hq = 0.1 kpc^Myr“^, r/p^max = 0.1. The shape of the background model 
can strongly influence the total lepton spectrum. 


For completeness, we wanted to test the synthesis model prediction against that 
obtained using detected radio pulsars, to ensure that the first is indeed higher, since it 


®http://galprop.stanford.edu/webrun/ 
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encapsulates both detected and undetected pulsars. We se 


F > 0, and d < 2 kpc from the ATNF Pulsar Catalog flManchester et al 


ected all pulsars with P < 0.1 s, 


200511 . We 


removed globular cluster pulsars, young pulsars (such as Vela), and known BW and RB 


systems by hanc 


(Shkloyskii 


, leaving us with ~ 80 MSPs. We used Shklovskii-corrected values for P 


1970) when available. We then repeated the calculation above, and plotted the 


result (not shown) in order to compare with that from the population synthesis (where we 
have ~ 100 sources within 1 kpc). We conhrmed that the “ATNF component” was lower 
than the “synthesis component”, as expected, since the detected pulsars should be a lower 
limit to the total number of sources predicted by the synthesis model. However, this is 
strongly dependent on e, with the “ATNF component” becoming closer to the “synthesis 
component” for higher e. This reflects the facts that the main contribution comes from 
nearby, powerful MSPs, and that the dominant contribution comes from pairs, the level of 
which very sensitively depends on pair multiplicity. 


Figure [16] shows the measured PF (e.g.. 


Accardo et al. 


l2014j) as well as the GALPROP 
and synthesis plus BW / RB contributions, for kq = 0.1 kpc^ Myr“^ and r7p,max = 0.1. 
The largest contribution is found ~ 100 GeV in the case of £ = 0.6 and F = 3.6 /iG. 
Figure [T7I is the same, but for hq = 0.1 kpc^Myr”^, r^p^max = 0.3, while in Figure [18] we use 
Ko = 0.01 kpc^Myr”^ ?7p,max = 0.1, and in Figure [IH] kq = 0.01 kpc^Myr“\ 77p,max = 0.3. 
In Figure [20l we show results for different choices of aj), the highest ratio (above 1 TeV) 
occurring for the lowest value of ctd- Figure [21] is for the background model of DIO, for 
Ko = 0.1 kpc^Myr”^, r^p^max = 0.1, while Figure [22] is for the background model of DIO, 
for Ko = 0.01 kpc^Myr”^, rjp^ranx = 0.3. We note that the BW / RB component makes a 
signihcant contribution at a few hundred GeV, increasing with r^p^max and decreasing with 
Ko, while the result is very sensitive to the choice of background model. Some parameter 
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combinations are excluded by the data, e.g., kq = 0.01 kpc^ Myr = 0.3, and 

£ = 0.6, depending on the choice of background model. 


Discussion 


Our results have shown that for certain ranges of parameters, MSPs could make 
a signihcant contribution to the local cosmic-ray lepton spectrum. On the other hand, 
our calculations show that some parameter combinations can also be ruled out. Shock- 
accelerated positrons and electrons from BWs and RBs make a much larger contribution 
than the rest of the MSP population (above ~ 100 GeV), since a much higher fraction 
of the pulsars’ spin-down power goes into particle power in these sources (compared to 
the pair cascades). The BW and RB contribution nearly reaches the observed positron 
fraction for hq = 0.01 kpc^Myr”^, T^p^max = 0.1, e = 0.6 and ReS = 3.6/iG. However, the 
extreme parameter combination kq = 0.01 kpc^Myr”^, r/p^max = 0.3 and £ = 0.6 can be 
ruled out, while kq = 0.01 kpc^Myr”^, -/^p^max = 0.1 and £ = 0.2 would predict a rise in 
both the electron and positron spectra and the positron fraction above 500 GeV. We also 
note that the uncertainty in the energy dependence of the spatial diffusion coefficient (and 
background model) may lead to large uncertainties in the flux predictions, especially at very 
high energies. Future measurements by AMS-02 extending the spectra to larger energies 
may help constrain even more of parameter space. We also note that the level of the pair 
spectra starts to saturate for larger e. This is because more photons convert to pairs at 
lower energies in this case (mostly due to the lower curvature radius of the distorted held 
lines). One therefore cannot increase the source hux without bounds by increasing e, since 
there is a limit to the maximum contribution one would obtain as a function of e. We have 
chosen values for e that are reasonabl e, simulating the range of of fsets obtained in newer 


solutions to realistic magnetospheres flHarding fc Muslimov 


2nilbll . 
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We found that the PF could be nicely reproduced if the population synthesis component 
were shifted to an energy higher by a factor of ~ 4 and at a level higher by a factor of 
~ 10. Uncertainties in the transport and source properties may account for some of this 
shortfall, so that these numbers are not too large. However, we rather interpret this as 
pointing to the fact that the synthesis component has the correct spect ral shape to explain 


Gaggero et ah 


20141) . given the assumed 


the data (i.e., hard enough spectral index; see, e.g., 
background model. 

We have taken the approach of computing the positron contribution from MSPs using 
reasonable parameter values rather than tuning the parameters just to £t the PF. This 
is because we do not expect the MSPs to explain all of the data, since there are many 
good arguments why young pulsars in PWNe, as well as SNRs, may make even larger 
contributions. The rise in PF may plausibly be the result of contributions by many sources 


(e.g., 


DiMauroetah 


20141) . One should therefore be careful not to overproduce the data 


by not considering the cumulative contribution from all viable sources. The argument may 
actually be reversed: if one has access to solid predictions for the contribution of several 
sources, one could in principle constrain parameters such as the shock acceleration efficiency 
so as not to overproduce the observed flux. In practice, this may be difficult, though, given 
the large number of free parameters and model uncertainties. 

It is clear that the predicted spectrum of secondary electrons and positrons from 
cosmic ray interactions in the ISM has a significant influence on 
PF, and therefore the need for a primary positron contribution. 


the rise anc 

shape 

Blum et ah 

(2OI3I 


that the PF fraction rise can even be explained exclusively by secondary positrons. Such a 
hypothesis would imply strong constraints on the properties of the primary positron sources, 
including the roughly 50,000 Galactic MSPs, sho ck acceleration in BW and RB binaries. 


and young PWNe. However, the calculations of 


Blum et af. 


(120131) have a factor 
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uncertainty, implying that there may be a primary contribution of equal strength around 
several hundred GeV, and even at a higher relative level at lower energies. Furthermore, a 
source that cuts off at lower e nergies will no t viola te this upper limit. We have also noted 


in Section [T] the conclusion by 


Dado fc Dail (120151) that a primary source of positrons is 


strongly required to explain the cosmic-ray data. We therefore envision room for both a 
primary and secondary contribution to the observed excess. 


It became clear during our study that the results are very dependent on the properties 
of the closest few sources, since they dominate any contribution from the larger, more 
distant, population. A change in the properties of these nearby sources for a new realization 
of the synthesis population may therefore impact our predictions. To investigate this 
matter, we split the population into two parts, and compared the cumulative electron and 
positron spectra from these two subpopulations. While there were some differences, the 
effect was fortunately rather minor, and we can therefore have some conhdence in our 
predictions. 


Our population synthesis uses radio survey sensitivity and Fermi three-year point 


source sensitivity maps, normalizing to the number of detected radio MSPs 
surve ys only, and to the detected 7 -ray MSPs in the Second Pulsar Catalog 


rom those 


Abdoetah 


20131 ). all of which are radio-loud. The MSPs discovered in radio followup observations of 


unidentihed Fermi sources are included in our simulated population of MSPs not detected 
by radio surveys but detected by Fermi as point sources. However, it is possible that there 
is a contribution to the cosmic-ray flux from BWs and RBs very close to Earth that, due to 
difficulty of radio detection in these eclipsing and/or obscured systems, have not yet been 
identified. Future detections will impact the normalization of the output from the synthesis 
component, and may also enhance the flux prediction from the BW / RB component. Our 
predicted contribution from the latter should therefore be seen as a lower limit, since we 
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have only considered detected BWs and RBs in our calculation. 


We note that shock acceleration in 


jinary systems s. 


that is modulated at the orbital period (lArons fc Tavani 


ionic 


1993 


ead to non-thermal emission 


Bogdanov et ah 


2011 1. None 


has been seen in the y-ray band so far, although there have been some det ections of X-r ay 


emiss ion modulated at the orbital pe riods in the BW sy stems B1957+20 flHuang et ah 


2012 1. J2215+5135, and J2256-1024 f Gentile et ah 2014). Future detections of such 
high-energy modulated signals will provide further conhrmation and constraints on the 
shock acceleration scenario. 


Kisaka fc Kawanakal (1201211 suggested that pair-starved MSPs may be responsible for 
a large peak in the total electron spectrum at 10 — 100 TeV, and that non-pair-starved 
MSPs with multiplicities of ~ 2 000 may contribute significantly (near 100%) to the PF 
above 10 GeV. Ther e are, however, a number of differences in our respective approaches. 
Kisaka fc Kawanakal (1201211 used hxed values for P and Bg for all members in their 
population. They furthermore assumed energy equipartition between the particles and the 
magnetic field, which seems to imply a conversion efficiency (from spin-down luminosity to 
particle power) of p ~ 50%, while we find rj ~ 1% from our pair cascade modeling. They 
also assume a lower average Galactic magnetic field {B = 1 pG). Finally, they integrate the 
injected spectra over the age of the MSPs while we follow a steady-state approach. Most if 
not all of these differences should lead to an enhanced particle flux in their case. 


It has been suggested that isotropy may be a discriminator between an astrou 


and dark matter origin of the rise in PF (e.g.. 


Biisching et al. 


2008a : 


Aguilar et al 


lysica l 


2014h . 


However, if s everal nearby source s contribute, any po 


out (see also 


Feng fc Zhang 


201511 . 


Linden fc Profnmo 


ential anisotropy may be washed 


(120131) also raised some issues, e.g.. 


inhomogeneous magnetic helds, d iffusion properties b eyond the standard assumptions. 


or pulsar proper motion. Indeed, 


Kistler et ah 


(120121) noted that turbulence in the local 
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magnetic field may lead to filamentary structures or “streams” of cosmic rays. These streams 
may redirect or wash out signatures of local sources that may otherwise have contributed 
under the assumption of isotropy (i.e., making the terrestrial spectra nearly featureless), or 
conversely lead to an enhancement in the contribution by otherwise negligible sources by 
concentrating their fluxes. One should lastly consider the “coherence length” or mean free 
path beyond which any anisotropy would disappear. Given these uncertainties, we would 
argue that anisotropy measurements may not be such a clear discriminator after all. 


7. Conclusion 

In this paper, we carefully assessed the contribution of MSPs to the cosmic-ray lepton 
spectra at Earth using a population synthesis code and a pair cascade code to calculate 
realistic source spectra. We also considered the contribution of binary BW / RB systems, 
which may further accelerate pairs escaping from the MSP magnetospheres in intrabinary 
shocks. 

We hnd that the predicted MSP particle flux increases for non-zero magnetic held 
offset parameters e. This is expected, since a larger value for the offset of the surface 
magnetic held with respect to the non-perturbed magnetic axis leads to an increase in the 
acceleration potential for some regions in azimuthal phase (and a decrease in others). This 
in turn results in an enhancement in both the number of particles (since the multiplicity 
will be higher) as well as the maximum particle source energy (given a larger local electric 
held in some regions). We hnd that the MSPs from the synthesis model make only a 
modest contribution to the terrestrial cosmic-ray hux at a few tens of GeV, after which this 
spectral component cuts oh. This is because the maximal injected particle energy is limited 
by a maximal electric held, which depends on the MSP source properties such as P, Bg, 
and e. The ehect of diherent Galactic magnetic helds and soft-photon energy densities is 



also shown: an increased field and densities lead to increased energy losses, and vice versa. 
We have bracketed these losses, and note that the fiux uncertainty is not too large, given 
these uncertainties in magnetic field and energy densities. Although the PF is somewhat 
enhanced above ~ 10 GeV, our added MSP synthesis component fails to reproduce the 
high-energy rise for the parameters considered. 

The BW / RB component contributes more substantially above several hundred 
GeV, given the fact that they further accelerate the electron-positron pairs in their strong 
intrabinary shocks. For some parameter combinations, this component may even exceed 
the measured positron spectrum, and may violate the PF at high energies, depending on 
the background model. 

Alternative sources of primary positrons such as young, nearby pulsars or SNRs should 
also contribute to the cosmic-ray electron and positron fiux. Future observations and 
modeling should continue to constrain the properties of these source classes, as well as 
improve our understanding of Galactic structure and particles within our Galaxy. 
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Fig. 1.— Histograms of period log,o(7^), period derivative log,o(7^), surface magnetic field 
log,o(Hs), and distance d characterizing the simulated present-day Galactic MSP population 
(Section 12.21) . 
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Fig. 2.— Comparison of several simula ted and measured pr operties of a population of 
detected radio-loud MSPs. Adapted from 


Gonthier et al. 


(120151). 
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Fig. 3.— Schematic of the fo r matio n of a shock upon collision of pulsar and companion 


winds. Adapted from 


Harding! (1990|) 
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Fig. 4.— Sample electron-positron pair spectra (number of pairs per second and energy) 
calculated for different periods P and offset parameters e, as indicated in the legend, and for 
a fixe d = 20 (i.e., = 2 ^ 10^ G). The x-axis indicates source energy in units of rrieC^. 


From 


Harding fc MuslimovI (2011b). 
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logio(E^^,/erg) log^g(E^^,/erg) 


Fig. 5.— Contour plot of log^o (F^min/-Fcut) vs. logio(M+) and logio(.Fcut) assuming P^as = 3 
and Bs = 5. Panel (a) is for ?7p,max = 0.1, and panel (b) is for ?7p,max = 0.3. Values near 
unity (dark red, i.e., the lower left corners) indicate that no solution could be found. 
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Fig. 6.— Comparison of transport done for the Galactic synthesis component (for one 
particular realization of the MSP population using the synthesis code) in the Thomson 
(TL) and KN limits. We compare Disk and Halo scenarios (which have the same value of 
Heff = 7pG in the case of the Thomson limit), and also consider results for different values 
of e and kq, as indicated in the legend. 
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Fig. 7.— Comparison of contribution of synthesis vs. BW / RB component, assuming 
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Fig. 8.— Same as Figure [71 but assuming kq = 0.1 kpc^ Myr“^ and ?7p,max = 0.3. 
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Fig. 9.— Same as Figure [71 but for i?efr = 3.6/rG, Uopt = 0.4 eVcm r^p^max = 0.1, 

kq = 0.1 kpc^ Myr“^, and ctd = 0.3 and 0.6. 
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Fig. 10.— Total MSP contribution (assumed to be equal numbers of positrons and electrons) 
to the leptonic cosmic-ray spectrum at Earth, assuming kq = 0.1 kpc^Myr”^ and r^p^max = 
0.1. Electron spectra appear at the top, while positron spectra appear lower down. The 
contribution from the synthesis component to the positron spectrum is visible at ~ 30 GeV 
(for £ = 0.6), and that of the BWs and RBs at ~ 1 TeV. The cool colors (purple, blue, and 


cyan) indicate spectra for £ = 0.0, 0.2, and 0.6. Green indicates th e “background” (non-Iv 


electrons and positrons using output from the GALPROP code 
standard parameters. Also show n are data from PAMEL A (red; 


(orange; 


Ackerman n et ah 


2012 1. and AMS— 02 (yellow; 


the website http://lpsc.in2p3.fr/cosmic-rays-db flManrin et al 


Vladimirov et ah 


Adriani et al. 


Aguilar et al. 


201111 for 


SP) 


2013 1. Fermi 


20141) . accessed via 


20i4 . 
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Fig. 11.— Same as Figure ITOl but for kq = 0.1 kpc^ Myr“^ and T^p^max = 0.3. 
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Fig. 12. 


Same as Figure ITOl but for kq = 0.01 kpc^ Myr ^ and ?7p,max = 
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Fig. 13.— Same as Figure ITOl but for kq = 0.01 kpc^ Myr ^ and ?7p,max = 0.3. 
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Fig. 14.— Same as Figure fTOl but for = 3.6/xG, Gopt = 0.4 eVcm“^, r^p max = 0.1, 
kq = 0.1 kpc^ Myr“^, and = 0.3 and 0.6. 
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Fig. 15.— Same as Figure ITOl but for the background model of DIO. 
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Fig. 16.— Measured fjAckerma.uu et a1. 


2012 


Adriani et al. 


2013 


Accardo et al. 


2014Ji and 


predicted PF (including “background” contributions from GALPROP in green and the syn¬ 
thesis plus BW / RB contributions from this work in purple, blue, and cyan, indicating 
e = (0.0, 0.2, and 0.6). Here, Kq = 0.1 kpc^ Myr“^ and r^p^max = 


0 . 1 . 
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Fig. 17.— Same as Figure [T6l but for kq = 0.1 kpc^ Myr ^ and r^p^max = 0.3. 
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Fig. 18.— Same as Figure [T6l but for kq = 0.01 kpc^ Myr ^ and ?7p,max 


0 . 1 . 
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Fig. 19.— Same as Figure [T6l but for kq = 0.01 kpc^ Myr ^ and ?7p,max 


0.3. 
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Fig. 20.— Same as Figure [T6| but for i?eff = 3.6 /iG, Gopt = 0.4 eVcm r^p^max = 0.1, 
kq = 0.1 kpc^ Myr“^, and = 0.3 and 0.6. 
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Fig. 21.— Same as Figure [161 but using the background model of DIO, for kq = 
0.1 kpc^ Myr“^ and = 0.1. 
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Fig. 22.— Same as Figure [161 but using the background model of DIO, for kq = 
0.01 kpc^ Myr“^ and ?7p,max = 0.3. 
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Table 1. Measured and derived parameters of BW pulsars. 


Name 

Pms 

P^ 

(10-20) 

(1004 ergs-1) 


d 

(kpc) 

Pb 

(h) 

-^^comp 

{Me) 

Oil 

^cut 

(TeV) 

Ref. 

J0023+0923‘= 

3.05 

1.15 

2.50 

4.88 

0.7 

3.3 

0.016 

1.01 

2.40 

1 

J0610-2100‘= 

3.86 

0.34 

0.36 

2.96 

3.5 

6.9 

0.025 

1.65 

3.04 

2 

J1124-3653‘= 

2.41 

0.57 

2.50 

3.05 

1.7 

5.4 

0.027 

1.40 

2.03 

1 

31301+0833'= 

1.84 

0.95 

9.36 

3.44 

0.7 

6.5 

0.024 

1.59 

1.37 

3 

31311-3430'= 

2.56 

2.08 

7.64 

6.01 

1.4 

1.56 

0.008 

0.61 

2.33 

4 

31446-4701'= 

2.19 

1.01 

5.93 

3.88 

1.5 

6.7 

0.019 

1.62 

1.52 

5 

31544+4937'= 

2.16 

0.31 

1.87 

2.12 

1.2 

2.8 

0.018 

0.91 

2.72 

6 

31731-1847 

2.34 

2.47 

11.9 

6.26 

2.5 

7.5 

0.04 

1.75 

1.23 

7 

31745+1017'= 

2.65 

0.23 

0.75 

2.02 

1.36 

17.5 

0.016 

3.07 

1.86 

8 

31810+1744'= 

1.66 

0.45 

6.08 

2.26 

2 

3.6 

0.044 

1.07 

1.86 

1 

31959+2048'= 

1.61 

0.72 

10.6 

2.80 

1.53 

9.2 

0.021 

2.00 

1.19 

9 

32047+1053'= 

4.29 

2.00 

1.56 

7.63 

2 

3 

0.035 

0.95 

2.78 

3 

32051-0827'= 

4.51 

1.23 

0.83 

6.14 

1 

2.4 

0.027 

0.82 

3.51 

2 

32214+3000'= 

3.12 

1.46 

2.96 

5.57 

1.32 

10 

0.014 

2.11 

1.59 

10, 11 

32234+0944'= 

3.63 

1.94 

2.50 

6.91 

1 

10 

0.015 

2.11 

1.66 

3, 5 

32241-5236'= 

2.19 

0.67 

3.90 

3.15 

0.5 

3.4 

0.012 

1.03 

2.12 

12 

32256-1024'= 

2.29 

1.58 

8.11 

4.96 

0.6 

5.1 

0.034 

1.35 

1.54 

1 


®'We have used the equation of state (EOS) described in Section 12.31 For canonical values, divide by a 
factor 1.56. 

'^We have used the EOS described in Section [231 For canonical values, multiply by a factor 0.78. 
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Fermi LAT pulsations have been seen from this pulsar. 

Note. — The columns are as follows: Pulsar name; pulsar period in milliseconds; intrinsic (Shklovskii- 
corrected) period derivative, as calculated from the announced spin-down luminosities; spin-down luminosity; 
surface magnetic field in units of 10® G; distance; binary period; companion mass; binary separation in units 
of 10^^ cm; spectral cutoff energy. We assume a pulsar radius of i? = 9.9 x 10® cm and moment of inertia of 
I = 1.56 X 10'^® gcm^. 


Ref eren ces. — (lilHessels et al 


(9) 


(120121): 151 


Keith et al. 


Fruchter et al 


( 20121 ) : 


(11990) : (10) 


(2011); (2) Espinoza et al. 


Bhattacharvva et al 


RobertsI ( 20111) : (11) 


Ransom 


( 20131); (7) 


Keith et al. 


(201 


3 ). 


( 2013); (3) Rav et al. ( 20121); (4) Pletsch et al 


( 120101 ) : ( 8 ) 


Barr et al 


(120131) : 
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Table 2. Measured and derived parameters of RB pulsars. 


Name 

Pms 

to 

0 

(1034 ergs-1) 


d 

(kpc) 

Pb 

(h) 

-^^comp 

(Mo) 

Oil 

^cut 

(TeV) 

Ref. 

J1023+0038 

1.69 

1.20 

15.4 

3.72 

0.6 

4.8 

0.2 

1.33 

1.33 

1 

J1628-3205 

3.21 

1.13 

2.11 

4.96 

1.2 

5 

0.16 

1.36 

2.15 

2 

J1723-2837 

1.86 

0.75 

7.18 

3.08 

0.75 

14.8 

0.4 

2.90 

1.09 

3,4 

31816+4510'= 

3.19 

4.03 

7.64 

9.34 

2.4 

8.7 

0.16 

1.97 

1.30 

5 

J2129-0429 

7.61 

43.54 

6.08 

47.4 

0.9 

15.2 

0.37 

2.94 

1.12 

6 

32215+5135'= 

2.61 

2.79 

9.67 

7.03 

3 

4.2 

0.22 

1.22 

1.55 

6 

32339-0533'= 

2.88 

1.39 

3.59 

5.21 

0.4 

4.6 

0.26 

1.30 

1.93 

7,8 


®'For canonical values, divide by a factor 1.56. 

'^For canonical values, multiply by a factor 0.78. 

Fermi LAT pulsations have been seen from this pulsar. 


Note. — The columns are the same as for Table [T] 


Ref erenc es. — (1) Archibald et al 


(1201,111 : (5) 


Kaplan et al 


(I2OI2I); ( 6 ) 


(200911: (2) 


Hessels et al. 


Rav et al.l ( 2012); (3) Robert: 


(j201l|); (7) 


Kong et al 


( 1201211 : ( 8 ) 


( 2011); (4) Crawford et al 


R.av et al 


(I 2 OI 4 I 1 






















































